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Abstract 

We consider a quantum system consisting of a one-dimensional chain of AI identical harmonic 
oscillators with natural frequency ui, coupled by means of springs. Such systems have been 
studied before, and appear in various models. In this paper, we approach the system as a 
Wigner Quantum System, not imposing the canonical commutation relations, but using instead 
weaker relations following from the compatibility of Hamilton's equations and the Heisenberg 
equations. In such a setting, the quantum system allows solutions in a finite-dimensional Hilbert 
space, with a discrete spectrum for all physical operators. We show that a class of solutions can 
be obtained using generators of the Lie superalgebra 0[(1|M). Then we study the properties 
and spectra of the physical operators in a class of unitary representations of g[(l|M). These 
properties are both interesting and intriguing. In particular, we can give a complete analysis 
of the eigenvalues of the Hamiltonian and of the position and momentum operators (including 
multiplicities). We also study probability distributions of position operators when the quantum 
system is in a stationary state, and the elfect of the position of one oscillator on the positions 
of the remaining oscillators in the chain. 
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1 Introduction 



In recent years quantum information theory has known an enormous expansion. This has boosted 
new interest in probabihstic and geometric aspects of state spaces of simple quantum systems. In 
this context, the dynamics of entanglement in a chain of coupled harmonic oscillators has been 
the subject of many papers One of the systems for which entanglement dynamics is 

being studied consists of a large chain of harmonic oscillators coupled by some nearest neighbour 
interaction In a popular model this coupling is represented by springs obeying Hooke's law. 
Then the Hamiltonian of the system is given by: 

-2 1 

k=l 

In other words, the quantum system consists of a string or chain of M identical harmonic oscillators, 
each having the same mass m and natural frequency uj. The position and momentum operator for 
the fcth oscillator are given by qk and pk', more precisely qk measures the displacement of the kth 
mass point with respect to its equilibrium position (see Figure 1). The last term in (jLljl represents 
the nearest neighbour coupling by means of "springs", with a coupling strength c (c > 0). Finally, 
we shall assume periodic boundary conditions, i.e. 

qM+i=qi- (1.2) 

Such quantum systems are also relevant in quantum optics (photonic crystals), or for describing 
phonons in a crystal E] . 
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Figure 1: Illustration of the quantum system: (a) the M masses in equilibrium position, (b) certain 
displacements q^ for each oscillator. 

In the standard approach for the quantum system governed by , one assumes the canonical 
commutation relations (OCR's) 

[qk,qi] = 0, [Pk,Pi] = 0, [qk,Pi] = i^hi- (1-3) 

Then, reformulating the problem in normal coordinates, the eigenstates of H can be described in 
some infinite-dimensional Fock space 

The approach of the present paper is more general. Instead of postulating the OCR's, we shall 
start from a more general quantization procedure. This procedure is based upon the compatibility 
of Hamilton's equations with the Heisenberg equations. Such systems are called Wigner Quantum 
Systems (WQS's) The idea is based on Wigner's observation - for the simple example of a 
one-dimensional harmonic oscillator - that this quantum system also allows solutions for which 
both Hamilton's equations and the Heisenberg equations are satisfied, but not the OCR's jHl- In 
other words, the OCR's are sufficient but not necessary conditions for Hamilton's equations and 
the Heisenberg equations to be compatible. Wigner's work led to the theory of parabosons and 
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parafermions in quantum field theory [Sl llUlITT] . and because of this attention its impact for ordinary 
quantum systems was somewhat overlooked. Another reason why WQS's did not receive immediate 
attention was because no general solutions could be constructed for the compatibility conditions 
of simple WQS's (apart from the one-dimensional harmonic oscillator). It was only much later 
- after the theory of Lie superalgebra's was completed - that Palev \12\ observed that classes of 
WQS-solutions for the n-dimensional harmonic oscillator are described by means of representations 
of the Lie superalgebras 05p(l|2n) and s[(l|n) or 0[(l|n). This algebraic or representation theoretic 
approach to quantum systems has revived the interest in WQS's |13 ( I14| [T5]. The WQS approach 
has so far been applied to simple systems of free harmonic oscillators, with some interesting and 
surprismg results [m ini CHI Cni EHI- Here it is, for the first time, apphed to a more realistic 
quantum system. 

In this paper, we shall study the system described by as a WQS. This implies that, apart 
from the standard solutions for which the OCR's hold, we shall also discover non-canonical solu- 
tions. In particular, we shall show that the quantum system allows solutions in a finite-dimensional 
Hilbert space, thus with a discrete spectrum not only for H but also for the position and momen- 
tum operators qk and pk- The class of solutions considered here is related to a g[(l|M) solution 
of the quantization or compatibility conditions, and we consider then a simple class of 3[(1|M) 
representations in which the physical properties of the quantum system are analysed. 

In Section [21 we analyse the compatibility conditions for the system (ll.lj) as a WQS. This leads 
to an expression of the Hamiltonian H in terms of a set of 2M operators (r = 1, . . . , M), which 
are themselves certain lineair combinations of the qk and pk- These operators should satisfy 
certain triple relations, see (|2.29|) . The task is then to construct operator solutions (constructions 
for the as operators acting in some Hilbert space) satisfying these triple relations. Although a 
complete set of solutions cannot be given, we show in Section |31 that certain generators of the Lie 
superalgebra 0[(1|M) satisfy the triple relations. In other words, we present algebraic solutions of 
the compatibility condition (or quantization condition) in terms of a Lie superalgebra. In order to 
have algebraic solutions that also satisfy the required unitarity conditions, the coupling constant c 
is bounded by a critical value cq, depending upon M. Since the problem has a g[(l|M) solution, the 
unitary representations of gl(l|M) serve as Hilbert space representations for the operators of the 
quantum system. In principle, all unitary representations of g[(l|M) are allowed for this purpose. 
In this paper, however, we shall concentrate on a particular class of unitary representations W{p), 
mainly for computational purposes. For these representations W{p), the actions of the 0[(1|M) 
generators, and in particular of the operators a^, are very simple expressions, see Section^ On 
the other hand, the class of representations W{p) is already sufficiently rich to exhibit intriguing 
physical properties of the quantum system under consideration. These properties are examined 
in the remaining sections. In Section 13 we study the energy spectrum of the quantum system in 
W{p). If c = (absence of coupling), the system consist of M independent (or free) identical one- 
dimensional oscillators, and it is easy to see that there are M + 1 equidistant energy levels in W{p), 
with multiplicities according to the level k {k = 0,1, . . . ,M). For < c < cq, these energy 
levels split: we can give a closed formula for the energy levels themselves, for their multiplicities, 
and for the total number of energy levels (which grows like 3^/^). Section El is devoted to the 
investigation of spatial properties of the chain of coupled oscillators, if the system is in one of 
the representations W{p). Olearly, since the representations considered here are finite-dimensional, 
the spectrum of position operators q^ and momentum operators pk is discrete. We manage to 
give the spectrum of these operators in closed form, but so far an analytic expression for a set of 
orthonormal eigenvectors is missing. Finally, we examine numerically for a simple example {M = 4) 
some position probability distributions of the oscillator system. The case of atypical representations 
W{p) (i.e. p < M — 1) \s examined in Section 7, followed by some concluding remarks in Section 8. 
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2 The quantization procedure 

In our approach we shall require that Hamilton's equations 



dH dH 

4 = ^, ^, = -^ (^ = 1,2,. ..,M) (2.1) 
and the Heisenberg equations 

i>k = \[HM, k = \[H,qk\ (A: = 1,2,...,M) (2.2) 

should be identical as operator equations. Since Hamilton's equations for the Hamiltonian 
take the explicit form 

Qk = —Pk, (2.3) 
m 

Pk = cmqk^i - m{u? + 2c) + cm^^+i, (2.4) 
the compatibility conditions read 

[HAk] = --Pk, (2.5) 
m 

[H ,Pk\ = —ihcm qk_i + ih'm{uP' + 2c) q^ — ihcmqk+i, (2-6) 
where k = 1,2, . . . , M, and - extending 1)1.21) - % stands for qM, or more generally 

Qk = Q kraodM, Pk=PkmodM (2-7) 

whenever k is out of the range {1, 2, . . . , M}. In other words, the task is to find operator solutions 
for gfc and pk such that the compatibility conditions H2.5() - H2.6|) . together with ()1.1() . are satisfied. 
Furthermore, since q^ and p^ correspond to physical observables, the operators should be unitary: 

Ql = Qk, p\=Pk {k = l,2,...,M). (2.8) 

Just as in the canonical treatment of the problem O E] j it will be useful to introduce "normal 
coordinates" and reformulate the problem in terms of these new coordinates. So let us consider 
the finite Fourier transforms of qt and pk'- 



Qr = ^^e- — gfc, (2.9) 

^ k=l 
M 

^ 1 \ — 2-Kirk , , 



The inverse relations are given by: 



qk 



Pk 



1 ^ 

^J^e^Q., (2.11) 



r=l 
M 



— y^e- — Pr. 2.12 



r=l 
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As for gfc and p^, see (|2.7|) . it will sometimes be useful to extend the indices by 

Qk = Q fcmodM, Pk = PkmodM- (2.13) 

Note that the unitarity conditions (|2.8|) imply: 

Ql = QM-r, Pl = pM-r (r = 1, 2, . . . , M), (2.14) 

and in particular, following the convention l|2.13() . Q\j = Qm and Pj.^ = Pm- An essential 
part is now to substitute (|2.11|) - (|2.12j) in the Hamiltonian H, given by (|l.lj) . and simplify this 
expression without assuming any commutation relations between the operators. For the term 
"YlkPk ~ YlkPkP^k^ one can, after the substitution, use the identity 

M 

E27Tik(r — s) 
e A? = M6rs, 

k=l 

and thus YlkPk ~ PrP) ■ Similarly, q"^ = QrQr- The coupling terms can be written as 

J2(^k - qk+i){qk - Qk+i)'^ = ^qkql + - ^{qkql+i + qk+iql)- (2.15) 

k k k k 

In the right hand side, the first two sums are of the same form as before; the last sum yields 

Y^{qkqi+i + qk+iql) = ^^{e'^QrQl + e^QvQl) 

k r 

= E 2^0^(^)441 (2.16) 

r 

Thus we obtain, just as in the canonical case lOj 

^ 1 2 

H = Y.[^Prn + '^QrQ^), (2.17) 
r=l 

where, for r = 1, 2, . . . , M, the quantities uJr are positive numbers with 

ojI = oj"^ + 2c - 2ccos(— ) = u;^ + 4csin2(— ), (2.18) 

and clearly 

^M~r = "-^r- (2-19) 

Substituting (|2.11l) - ()2.12|) in H2.5|) - (|2.6() yields the set of compatibility conditions for the new oper- 
ators: 

[H, Qr\ = --PI [H, Ql] = --Pr, (2.20) 
m m 

[H, Pr] = ihrriLolQl, [H, PJ] = ihmujf.Qr. (2.21) 

The task is now reduced to finding operator solutions for Qr and Pr such that the compatibility 
conditions together with (|TT7|) . are satisfied. 
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As a final step it is convenient to introduce linear combinations of the unknown operators 
and Pr, say and a~ (r = 1, 2, . . . , M), by 



+ ' py (2.22) 

2h ^j2mL0rfi ^ 



mujr 



2h ^/2moJrfi ^ ^ 

with 

(a^)t = a^. (2.24) 
Observe that the inverse relations take the form 



Qr = ^^yM-r + ar), r = 1, . . . , M - I; QM = ^^{ai, + al,); (2.25) 

Pr = i^l^^{at-al,_^), r = l,...,M-l; = i^^^{a+ - a]^,); (2.26) 

with similar expressions for qI and Pj . In terms of the new set of unknown operators (r = 
1,2, ... , M), the Hamiltonian (|2.17|) becomes: 

H = '^^{a-,a^} = J2^io,rat + a+a^). (2.27) 

r=l r=l 

It is essential - and the reader should verify this - that in going from 1)2. 17(1 to ()2.27|) . no com- 
mutation relations among the operators Qr and Pr are used, but only identities like (|2.1l'{j) . (|2.14|1 
and (HHH). A final and simple calculation, using ((T^ and shows that (jT^ - l^^ is 

equivalent to 

[H,af] = ±hiOraf, (r = 1,2... ,M). (2.28) 

Thus 

Proposition 1 In the approach of as a Wigner Quantum System, the problem is reduced to 
finding 2M operators af- (r = 1, . . . , M), acting in some Hilbert space, such that (a^)^ = and 

M 

[J2^j{ajay +ayaj),at]=±2L0raf, {r = 1,2 . . . , M). (2.29) 

i=i 

The operators corresponding to physical observables and pk are then known linear combinations 
of , and the Hamiltonian H is given by \2.21\) . 

As we shall see in the following section, this is an algebraic problem that has a class of solutions in 
terms of the Lie superalgebra 0l(l|M). 

Before concentrating on an algebraic solution, let us end this section with a few words about the 
time dependency of the operators. The time dependency of qt and pk is determined by (|2.H|) - (|2.4j) . 
From these equations and (|2.9() - H2.1U|) it follows that 

Qr = — Pr = -mujlQl. (2.30) 
m 

Using ^TF^-^T^ yields 

a~ = —iiOrCLr , CL^ = iiOrtt'^ , (2.31) 
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with the evident solution 

a±(t) = e±''^'-V(0). (2.32) 

So it is sufficient to have solutions for the operators at time 0, = a^(0). The time dependence 
for a^(t) is given by (|2.;-i2j) . and since all operators of the quantum system can be expressed in 
terms of a^, their time dependence follows. As a consequence, we shall concentrate on solutions of 
the system (fT77|l - (fT^ at time t = 0. 

For completeness, it should also be mentioned that if the OCR's (|1.3|) hold, then the operators 
satisfy the usual boson relations [a^,a^] = 0, [a~,a^] = Srs- In that case (|2.28|) follows 
automatically from H2.27() . 



3 Algebraic solutions of the compatibility conditions 

The set of relations (|2.29|) . together with the conditions (|2.24j) . are reminiscent of the algebraic 
relations satisfied by a set of 0l(l|M) generators We shall show that our problem has indeed 
a solution in terms of the Lie superalgebra g[(l|M) or s[(l|M). Let us first recall the definition of 
g[(l|M): it is a Lie superalgebra with basis elements Cjk, with j,k = 0,1, ... , M. The elements Cko 
and eok {k = 1, . . . , M) are odd elements, having degree deg(efco) = deg(eofc) = 1; the remaining 
basis elements are even elements, having degree 0. The Lie superalgebra bracket is determined 

by [HI Ea EH] 

In a representation, or in the enveloping algebra of £|[(1|M), the bracket (where x and y 

are homogeneous elements of 0[(1|M)) stands for an anti-commutator if x and y are both odd 
elements, and for a commutator otherwise. The Lie superalgebra s[(l|M) is the subalgebra of 
g[(l|M) consisting of elements with supertrace 0, or also s((l|M) = |0l(l|M), g[(l|M)]. 

For a Lie superalgebra one can also fix a star condition, i.e. an anti-linear anti-involution. For 
g[(l|M) or s[(l|M) such a star condition is fixed by a signature a, i.e. a sequence of plus or minus 
signs a = (cJi, . . . , aM) and 

(eok)^ = akCko, {k = l,...,M) (3.2) 

thus with each equal to +1 or —1. We are particularly interested in the case where all a^s are +1 
since this corresponds to the "compact form" u(l|M) of g[(l|M) |24j . for which finite-dimensional 
unitary representations exist j25j . 

We shall now show that there exist solutions of the form 

2~ / 2 

afcCfco, at = \ — ctl^keok, {k = l,...,M) (3.3) 

LOk V t^fc 

with Qfc certain complex constants to be determined. First of all, note that by (|3.2() . the unitarity 
condition ()2.24|) is automatically satisfied. With 1)3. 3() . the Hamiltonian 1)2. 27() becomes 

M M 

H = hi^'^ajlajl'^^ eoo + /i^o-fc|afcP Cfefc, (3.4) 
j=i k=i 

and the commutator of the above with (|3.3|) yields: 

M 

[iJ,a^] = ±;i(^cTj|aj|2-cTfc|afc|2)a± {k = l,...,M). (3.5) 
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These should coincide with the compatibility conditions 1)2. 28|) . Thus we get a system of M equa- 
tions in the unknown coefficients a^: 



M 

^aj\aj\^ - ak\ak\^ = ujk (A: = 1, 2, . . . , M). (3.6) 
It is easy to verify that a solution for this set of equations is determined by 

M 



—^Y.^r (3.7) 



For further use, it will be convenient to introduce the following numbers: 

M 



with LOj the values fixed by (|2.18j) . Note that 

M M 



(3M-k = Pk, P='^Pj='^^j^ (3-9) 



and thus H can be rewritten as 

M 

H = h{Peoo + ^l3kekk)- (3.10) 

k=l 

Remember that we are primarily interested in the signature with all cjfc's equal to +1. Since 
Cfclc^fcP = the question is whether there exist solutions such that all /3fc's are positive. At first 
sight, (|H.8j) indicates that (3k is equal to —Uk plus "some average value" of the Wj's, and hence 
one would expect half of the P^s to be negative and half of them to be positive. We shall show, 
however, that under certain conditions ("weak coupling", i.e. a small value for c), all /3fc's are indeed 
positive. First of all, note that for c > 0, 

Pi > (32> ■■■ > /3[A//2J ) PlM/2\ < PlM/2\+l < ■ ■ < Pm, (3-11) 

since a similar property holds for the values ujk- Thus if /?[m/2J > 0) then all P^s are positive. The 
value of PiM/2\ depends on the value of c; if c = then indeed all Pk = uj/{M — 1) are positive. 
So by continuity as a function of c there will be a certain interval ]0, co[ where all P^s are positive. 
This critical value cq is the c- value for which P[m/2\ = 0- For general M, this is a complicated 
transcendental equation that can be solved only numerically. Table 1 gives the numerical solutions 
for this equation, for M ranging from 4 to 21 (for M = 2 or M = 3 the PkS are always positive). 
The following proposition gives a lower bound for the critical value cq, such that all /S^'s are positive. 

Proposition 2 An upper bound for c is determined by: 

< c < ^j^zr^ ^ Pk>0, forl<k<M. (3.12) 
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Table 1: Critical values cq/uP 



M 


u / 


M 


u / 


4 


0.9873724357 


13 


0.10546881460 


5 


0.7500000000 


14 


0.09256321610 


6 


0.3457442295 


15 


0.08687882025 


7 


0.2982653656 


16 


0.07814800074 


8 


0.2061705212 


17 


0.07388896853 


9 


0.1851128402 


18 


0.06760983697 


10 


0.1464642846 


19 


0.06429500840 


11 


0.1343028683 


20 


0.05957194222 


12 


0.1134651313 


21 


0.05691629341 



Proof. By definition (|3.8j) , /3fc > if and only if 



1 ^ 



M 



[M - ifuol < Y,^j + 2 J^WiWj. 

j=l i>j 



(3.13) 



We write the right hand side of (|3.13|1 as a series with respect to c. For the first sum, we find: 
E^l = E ( + 4^^-2(5) ) =Mu;' + 2cY^(l- cos(^) ) = M{u;' + 2c). 



j=l j=l V 

Here, we used the Lagrange identity: 



M 



1 ^ sin((Hli+i)£) 



1 + cosfx) + cos(2x) H h cos(nx) = - + „ . 

2 2sin(|) 



In order to evaluate the second sum in (|3.13() . let 



i>j 



i>j 



F{c)=Y^ ^i^i = E \P + 4csin2(^)-^w2 + 4csin2(:^). 



It is then clear that 



F(0) = to' 



;M(M- 1) 



One can write the derivate of ujk with respect to c as: 



dwfc 2 o/k-K 



and one finds that 



since a;^ > 0. This implies that F{0) < F(c) for c > 0. Thus it follows from (|3.13|) that 



M 



{M-lfujl<J2^] + 2F{0) 



(3.14) 



(3.15) 



is a sufficient condition for (3^ > 0. Since w| < w^^y2j < + 4c for 1 < /c < M, it is sufficient to 
solve the following inequality 



(M - lyiuj^ + 4c) < Af(a;^ + 2c) + uj^M{M - 1) c < 



a; 



2 



2(M - 2)^ 



leading to (IrTT^ . □ 
Note that for M = 2 or M = 3 there are no conditions: for M = 2, /3i = 0^2 > and 

/32 = > 0; for M = 3, /?i = /32 = w/2 > and /?3 = (w^ + 3c)i/2 _ > 0. 
We can now summarize the main result of this section in the following 

Proposition 3 For fixed M , let c satisfy 

(no condition if M = 2 or M = 3). Then the compatibility conditions \2.27\ )- ll2.28\) have a solution 
for the operators in terms of q[{1\M) generators: 



o-k = \— ^ko, = 1/ — eofc, {k = l,...,M) (3.16) 

with (3k given by h3. t^) . The unitarity conditions i2.24\ ) are equivalent with the star condition 

(eofc)''" = Cfco- (3.17) 



4 A class of 0[(1|M) representations 

In order to study properties of the given quantum system related to the 0l(l|M) solution of the 
previous section, one should consider representations of 0[(1|M) for which 1)3. 17() holds. These are 
known as the unitary representations (or star representations), and have been classified by Gould 
and Zhang [2^. For the explicit actions of the g[(l|M) generators on a Gel'fand-Zetlin basis for 
these unitary representations, see These actions becomes fairly complicated, so in this paper 
we will concentrate on a particular class of representations, the so-called Fock type representations 
Wip). 

Without going into the details of the construction of such representations, we briefly summarize 
their main features here. Further details can be found in |27[ I28j 

The representations W{p) are labelled by a number p, with either p G {0, 1, 2, . . . , M — 1} or 
else p G M with p > M — 1. We describe the representation by giving the basis vectors of the 
representation space W{p) and the action of the 0[(1|M) generators on these basis vectors. 

For p G {0, 1, 2, . . . , M — 1}, the basis vectors of W{p) are given by: 

wie) = w{6i,e2,...,eM), eie{o,i}, and\6\ = 6i + --- + eM <p- (4.i) 

Note that in this case the dimension of W{p) is given by 

dimW{p) = j2(^)- (4-2) 

fc=o ^ ^ 

For p real and p > M — 1, the basis vectors of W{p) are all vectors w{9) = w{6i, 62, ■ ■ ■ , 9m) 
with each 0i S {0, 1}. Clearly, for p > M — 1, the dimension of W{p) is given by 2*^. 
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The action of the 0[(1|M) generators on the basis vectors of W{p) is now given by: 

eoow{d) = {p - \e\) w{e); (4.3) 

ekkw{d) = ekw{e), {l<k<M)- (4.4) 

ekow{e) = (1 - ek){-if'^-+''"Wv-\o\ w{ei, . . . , 0,, + 1, . . . , Om), {i<k< m)- (4.5) 



e^kwiO) = ek{-ir^-^''-Wp-\G\ + ^w{eu. . . , - 1, . . . , Om). {l<k<M). (4.6) 

The action of the remaining 0[(1|M) basis elements can be determined from the above formulas, 
and one finds (for 1 < j < k < M): 

ejkw{d) = Bkil - ej^-lf^^-^'"-^ w{e^, ...,e, + l,...,dk-l,..., 9m), (4.7) 
ekjw{e) = -9j{l - 9k){-lp+-+'''-' w{9,,. . . , % - 1, . . . , ^fc + 1, . . . , ^A/). (4.8) 

With respect to the inner product 

{wmH0'))=6e,e', (4.9) 
the representation W{p) is unitary for the star condition 

The representations W{p) with p > M — 1 are typical irreducible g[(l|M) representations; 
those with p £ {0, 1, 2, . . . , M — 1} are atypical irreducible g[(l|M) representations H^lUnj- In 
Sections 5 and 6 we shall develop results for the typical representations; the atypical case is treated 
in Section 7. 

Note, by (|3.16j) and (|4.5j) - (|4.()j) the action of the operators on the basis vectors w{9) of W{p) 
(fc = 1, . . . , M) is given by: 



a^w{9) = J ^(1 - 0,)(-l)^i+-+e'^-i ^^(^1, ...,9k + l,...,9M), (4.10) 

V 



«(^) = J^9k{-lf^+-^''-^ Vp-\0\ + 1 w{9i, ...,9,-1,...,9m). (4.11) 

So the operators raise or lower 9k by one unit (if allowed). This means that, for p > M — 1, the 
basis vectors w{9) of W{p) have a Fock basis construction, by letting |0) = wi^l, 1, . . . , 1). Then 
a-|0) = and 

w{9) ~ iat)'-'H4)'-'' ■ ■ ■ (4-12) 
This is the reason why W{p) is referred to as a Fock representation. 

5 On the spectrum of H in the representations W{p) 

For any p > M — 1, the representation W{p) is of dimension 2^^. Under the solution (|3.16|) . the 
Hamiltonian H takes the form p.lUf) . 

M 

H = h{f3 eoo + '^Pk e-kk)- 

k=l 
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with /3 = Ylh=i h = Ylk=i ^k- Since the actions of etk (fe = 0, 1, . . . , M) are diagonal in the basis 
w{9), see (|4.3p - (|4.4() . it follows that the vectors w{9) are eigenvectors for H: 



Hw{e) = hEew{9), (5.1) 



with eigenvalues 



^ M - 2 ^'^ 

Ee = Pip -\0\)+Yl = Pip - - E ^k^k- (5.2) 

k=l k=l 

In this expression, 9 = (^i, . . . , 9m), with each 9k G {0, 1}, and \9\ = X^fcli ^k- 

In the case without coupling (c = 0), all PkS are the same: Pk = jfzj P = m-i ^- The 
eigenvalues of H are then 

huj{p- 



M -1 



The multiplicity of this eigenvalue is i^ji)- In other words, there are M + 1 distinct energy levels, 

M 



equally spaced with steps of unit huj. The lowest energy level corresponds to = ^Pjj 

coM, and the highest to -E'(o,...,o) = ^Ptj~\_- 

We are mainly interested in the weak coupling case (c > 0, but c < cq). Also in this case, it is 
easy to describe the energy levels through (|5.2|) . but the analysis of their multiplicity requires some 
further attention. For this purpose, observe that by ()3.9() 



M M M-1 M 

(^kPk 

k . k=l k=0 k=l 



kPk — "^^QkPM-k — E (^M-kPk — '"^^M-kPk, (5.3) 



where in the last step we have followed the convention that Po = Pm, and we have set = ^Af- 
Let wi9) be an arbitrary eigenvector H with eigenvalue HEq. Obviously, by (|5.2|) and (|5.3j) all basis 
vectors 'wi9') which are obtained by swapping 6i and 9M-i for arbitrary indices i yield the same 
eigenvalue hEg. The multiplicity of this eigenvalue is thus (at least) 

Further inspection of the last expression in (|5.2|) shows that this is indeed the multiplicity of the 
eigenvalue, for < c < cq. 

In Figure 2 we give a plot of the energy levels for M = 4 and M = 5, as an illustration of the 
above. As c increases, the M + 1 equidistant energy levels for c = split up in different levels, with 
smaller degeneracies. 

It is worth noting that one can say something extra about the number of energy levels, for 
arbitrary M, also when c ^ 0. For M = 4, the five levels 1 + 1 + 1 + 1 + 1 at c = become 
1 + 3 + 4 + 3 + 1 levels for c > [just counting the energy levels, disregarding their multiplicities], 
see Figure 2(a). We list here the number of levels for a few M- values, when c > 0: 



M 


number of levels 


total number 


1 


1+1 


2 


2 


1+2+1 


4 


3 


1+2+2+1 


6 


4 


1+3+4+3+1 


12 


5 


1+3+5+5+3+1 


18 


6 


1+4+8+10+8+4+1 


36 


7 


1+4+9+13+13+9+4+1 


54 


8 


1+5+13+22+26+22+13+5+1 


108 
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(a) M = 4 



(b) M = 5 




Figure 2: (a) The energy levels of the quantum system for M = 4 in W{p), where we have taken 
h = uj = 1, p = M, and c varies from to cq. The vertical axis gives the energy values. The 
numbers next to the levels refer to the multiplicity. When c = there are M + 1 = 5 energy 
levels, with multiplicities (1,4,6,4,1). When < c < cq, the energy levels split up in 12 levels, 
with multiplicities 1 or 2. (b) The same illustration for M = 5. There are 6 levels for c = with 
multiplicities (1,5,10,10,5,1), and there are 18 levels for < c < cq, with multiplicities 1, 2 or 4. 

Let T{M,k) {k = 0,1,..., M) be the number of levels "per split" (the numbers in the middle 
column), and L[M) = '^kLoT{M, k) be the total number of energy levels for c > 0. T{M,k) is 
also the number of distinct energy levels for all 9 with |^| = k. Analysing this, using (|5.4|1 . one 
finds 

T(M, A;) = r(M -2,k)+ T{M - 2, A; - 1) + T(M -2,k-2). (5.5) 
Summing now over all k yields a simple recursion for the number of levels L{M): 

L{M) = 3L{M -2). (5.6) 

So we obtain the following result for the number of energy levels when < c < cq, depending on 
whether M is even or odd: 

L(2n + 1) = 2-3", L(2n + 2) = 4-3", (n = 0, 1, 2, . . .). (5.7) 

Note that the numbers L{M) belong to a known sequence, see entry A068911 in (HOI, with a simple 
generating function: 

G{x) = f; L(M) x*^ = ii^. (5.8) 
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6 On the spectrum of position operators in the representations 
W{p) and spatial properties 

The purpose of the present section is to consider some geometric aspects of the given quantum 
system. In first instance, we shall analyse the spectrum of the position operators qr in the repre- 
sentations W{p). By construction, these Hilbert spaces W{p) are finite dimensional, so the position 
operators have a discrete spectrum. 

The determination of the spectrum for general M and p requires a lot of work. This goes in 
two steps: first the operators will be considered; this is the difficult part. Once the spectrum of 
is analysed, that of qr follows rather easily. The spectrum of the momentum operators p^ and pr 
is completely similar. In fact, we shall see that the structure of and q^ as operators in 0[(1|M) 
is equivalent, and they can be treated simultaneously. 

Note that, due to the symmetry of the system and (|1.2j) . the spectrum of q^ and qr will 
be independent of r. 

6.1 Eigenvalues and eigenvectors for 

We start by writing the operator q^ in terms of the gl(l|M) basis elements. Using (|2.11|) . one finds 

Qr = ^{Qr,ql} = ^Z]2^e^^^^^'3fc,0^}• 
fc=l 1=1 

Using ()2.25|) and the solution ()3.16|) . this yields 

f- M MM I r, /J 

9r- = ^ttI — eoo+ > > e m e/fc). 6.1) 

k=l k=l 1=1 

In a similar way, one finds 

, M MM 

Pr = ^{E(3>^eoo + Y.E • (6-2) 

k=l k=l 1=1 

As operators of g[(l|M), these two expressions are structurally equivalent: (Jr is obtained from the 
expression of Pr by formally replacing every parameter /3k by /J^/w^, and multiplying by an overall 
factor l/rm?. So their spectral analysis is equivalent. Since the expression of is somewhat simpler 
from the point of view of notation (no denominators in the factors), we shall first concentrate on 
the analysis for p,^. 

Following (|6.2() . it will be useful to introduce the even 0((1|M) operators 

M M 

X — ^ \ — ^ 2'Kir{j — k) J 

Sr = 2^2^e M ^yPjPkejk■ (6.3) 
k=i j=i 

The main ingredient of our analysis is the following technical lemma: 

Lemma 4 In the representation W{p), the operator Sr has only two eigenvalues, namely (3 = 
X^fcli /5fc O'^d 0, each with multiplicity 2^^~^ . 

Proof. Consider a vector of the form 

M 

e— v^tx;(l'=), (6.4) 

k=l 
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where w{l^) denotes the basis vector of W{p) with |^| = 1 and the "1" occurring in position k 
(counting from left, and starting from one). For instance, if M = 4 then 

wil^) =w{l, 0,0,0), wil"^) = w{0, 1,0,0), = u;(0, 0,1,0), ^(1^) = ^(0, 0, 0, 1). 

We extend the notation to w{l^^ • • • l'^") which denotes the basis vector with \9\ = n and with "l"'s 

in the positions ki up to kn- 
it is easy to verify that Vr is an eigenvector of Sr with eigenvalue f3. 
In g[(l|M), one can see that 

M 

I 27virl ^ — ^ , 27Tirj 

[eio,Sr] = -VAe A/ 2^y/3je m ejo, (6.5) 
and thus the action 1)4. 5() implies that 

[eio,Sr]Vr = 0. (6.6) 

This means that e/o^r (provided it does not vanish) is also an eigenvector of Sr, with the same 
eigenvalue /3. Our goal is now to show that 

ez„oe«„_iO • • • ei.ovr (6.7) 

is also an eigenvector of Sr with the same eigenvalue f3, provided all li are different (and different 
from M). If one would apply the same e;o twice, the resulting vector vanishes since {eio, eko} = 0, 
by p. If) . This means that we can also assume, without loss of generality, that M > In > ln~i > 
■■■>li>0. 

Using ()4.5|) . one can write 

ei„oei„_,o---ei,oVr^Y.^-lY VKe'-^w{l''---l'nH''+'---l'"), (6.8) 

t=0 k=h+l 

with Iq = and In+i = M + 1. Using (|6.5|) . this implies that 

M n k+i-i 

ho, Sr]ei^oei„.,o ■ ■ ■ ei.oVr ~ Yl V^^t J^(-l)* ^ v^e^e,-ou;(l'i • • • l'n*^l''+i • • • l'"). 

j=l t=0 k=lt+l 

(6.9) 

Keeping in mind the action of Cjo on a basis vector w{9), one sees that this expression is a linear 
combination of basis vectors w{6) for which \9\ = n + 2, having a "1" in the positions h, I2 up to 
In and in two extra positions x and y. Consider such a vector 

y;^!^! . . . \hx \^\lrx + ^ . . . \hy 1^ I'iy + l • • • l'" ) . 

In expression (|6.9|) this vector will appear twice, once with k and j playing the role of x and y 
respectively, and one vice versa. In the first case the coefficient of this vector is 

^p-{n + l)^/M'ye^^^ X (-1)^^ X {-lYy+\ 
while in the second case it is 

Vp-(n + l)v^e^^^ X (-1)*^ X (-1)^^. 
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Since these two coefficients cancel, we can conclude that 



[e/o, Sr]e/„oQ„_iO • • • ehoVr = 0. (6.10) 

Following (|6.6|) and (|6.1()|) . all vectors of the form (|6.7j) are eigenvectors of Sr for the eigenvalue /3. 
It remains to find the number of linearly independent eigenvectors, i.e. the multiplicity of the 
eigenvalue (3. 

For n fixed, consider the (^^~^) vectors ()6.7|) with 1 < h < I2 < • • • < In ^ M — 1. Expressing 
these in the w{6) basis by (|6.8|) . in total basis vectors are involved. Each basis vector w{6) 

occurs {^^^) = n + 1 times, except those vectors that have a "1" in position M. There are 
C^n^) ^^^^ vectors. Next, let A be the {^~^) x matrix, consisting of the coefficients of 

the {^~^) vectors 1)6. 7p written in terms of the basis vectors. Select in A those columns 

that correspond to basis vector having a "1" in position M. This submatrix is equivalent with a 
diagonal matrix of which the diagonal elements are proportional to \/ (3m- So, the vectors (|6.7|) 
with l</i<---</„<M — 1 are linearly independent. 

Furthermore, it is immediately clear that the vectors 

ei„oei„_iO • • • ei^o^r and e^Qei^_^o ■ ■ ■ ei^oVr 

are linearly independent if n 7^ t (they are linear combinations of basis vectors with different |^|). 
The conclusion is that we have found 

Err)-"-' 



n=0 



linearly independent eigenvectors of Sr with eigenvalue /3. Note that for some fixed |^|, there are 
linearly independent eigenvectors of Sr with eigenvalue /?. 

The eigenvalue of Sr also has multiplicity 2*^~^. This is seen in a completely similar way, 
starting with the vector 

M 

k=l 

and acting repeatedly with eg/ (with 1 < / < M — 1) on this vector. Herein we have extended the 
notation of (|6.4j) : w{0^) denotes a basis vector where every 6j = 1 (j 7^ k) except dk = 0. For fixed 
1^1, there are thus linearly independent eigenvectors with eigenvalue 0. □ 

We can now describe the eigenvalues of p^. Since 



eoowie) = {p-\e\)w{0), 

it follows from (|6.2j) that an eigenvector of Sr which is a linear combination of basis vectors with 
fixed 1^1 will also be an eigenvector of p"^. The eigenvalues of p'^ are thus given by ^{p — K)P, 
for < K < M — 1, with multiplicities 2(^^^). More in particular, the eigenvectors of p^ with 
eigenvalue ^{p — K)j3 arise in two ways, one set having \6\ = K + 1 and containing vectors of the 
form 

The other set has \9\ = K and contains vectors of the form 
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Unfortunately, the vectors (for a given eigenvalue) constructed here are not orthogonal. For any 
fixed M, one can construct a set of orthogonal eigenvectors by Gram-Schmidt orthogonalisation. 
But so far we cannot give a closed analytic expression for some set of orthogonal eigenvectors. 

Let us also state the result for the squared position operators q^. Recall from the beginning of 
this subsection that every /3fe should be replaced by j3k/oj'j.-i thus (5 should be replaced by 

7 = E^- (6-11) 

Proposition 5 The operator has M distinct eigenvalues given by xj^ = ^^(p — K)'y, where 
< K < M — 1. The multiplicity of each x\ is 2(^'^^). The eigenvectors of q^ with eigenvalue 
xj^ arise in two ways: there are (^'^^) vectors with \9\ = K + 1 of the form 

and there are = (^j^"*^) vectors with \6\ = K of the form 

Herein, 

ur = Ve^^u;(l^), ur = Ve-^(-l)™t/;(0^). (6.12) 
6.2 Eigenvalues for qr 

We have shown that the eigenvectors of q'^ with eigenvalue xj^ have either |^| = i^or \9\ = K + 1. 
Let '4'r,XK be an eigenvector of qr with eigenvalue xk- Such an eigenvector necessarily has the form 

'>Pr,XK = X] ^d,r,XK'w{0) + E Ce,r,XKw{0), (6.13) 

\e\=K \e\=K+i 
with C some constants. Thus one can write 

Ce,r,XKQr'w{9)+ E Ce,r,XK^rW{e) = XK ^ Cq ,r,XKw{0) +X K ^ Ce ,r,XK'^{0) ■ (6.14) 

\e\=K \e\=K+i \d\=K |6'|=_ft:+i 

But the action of qr on a basis vector w{9) is necessarily a linear combination of basis vectors w{9') 
with \9'\ = 16*1-1 and \9'\ = \9\ + 1: this follows from (tTTTl) . (IT^ and KT^ - ^jTH^ . Thus the 
first sum on the left hand side of H6.14() is a linear combination of basis vectors with \9\ = K — 1 
and \9\ = K + 1, while the second sum on the left hand side of ()6.14|) is a linear combination of 
basis vectors with \9\ = K and \9\ = K + 2. Of these four linear combinations the first and the last 
vanish (since they do no occur on the right hand side), and it follows that 

C0^r,XKQrW{^) = ^ C'e,r,a;x^(^) • (6-15) 

\e\=K \e\=K+i 
Combining this with H6.14() implies 

Ce,r,XKQT'w{9) = XK Ce,r,XKw{0)- (6.16) 
\e\=K+i \e\=K 
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We will now show that, given H6.13() . 



V'r-xx = X] ^e,r,XK'>^{d) - X] Ce,r,XKwi0) (6.17) 

\9\=K \e\=K+i 

is an eigenvector of qr with eigenvalue —xk- In fact, the action of qr on ip^-XK follows directly 
from (|6.15|) and (|6.16|) and yields qri^r-XK — ~XK'4'r,-XK- 
Thus we have shown: 

Proposition 6 The operator q^. has 2M distinct eigenvalues given by ±xk = :;^{p ~ 

where < K < M — 1. The multiplicity of the eigenvalue ±xk is {^^x^)' '^^^ eigenvectors of qr for 
the eigenvalue ±xk contain, when expanded in the standard basis w{6), only vectors with \9\ = K 
or\9\ =K + l. 

So far, we have no simple analytic expression of the (orthogonal) eigenvectors of qr in terms of 
the standard basis vectors w{6). 

6.3 Position probability distributions for stationary states 'w{9) 

Consider the eigenvectors '4^r,x,g of the position operator qr for the eigenvalue x expanded in the 
w{6)-hasis: 

'>Pr,x,g = '^C0^r,x,gW{9), (6.18) 

e 

and assume that these vectors are orthonormal, i.e. 

{iJr,x,g\i'r,y,h) = ^Ce,r,x,gCe,r,y,h = ^x,y8g,h- (6.19) 

e 

In the above, g (or h) stands for a multiplicity label for vectors with the same eigenvalue x: e.g. 
when X = ±xk then g runs from 1 to (*^^). We know from the previous subsection that for 
X = ±xk, only coefficients Cg^r,x,g with 16*1 = i^'or \6\ = K + 1 will appear. 

Let us now suppose that the quantum system is in a fixed eigenstate w{9) of H (a stationary 
state). The expression 

P{9,r,±XK)= Yl \Co,r,±x^,g\^ (6.20) 

5=1 

yields the probability of "measuring" the value ±xk for the position of the rth oscillator when the 
system is in the state w{9). Plotting all these values P{9, r, ±xk) for K = 0, . . . , M — 1 yields the 
probability distribution of oscillator r in the stationary state w{0). 

It will be interesting to look at an explicit example of such probability distributions. First of 
all, due to the earlier mentioned symmetry of the system, these probability distributions will be 
independent of r; so we need to plot it for one r- value only (say r = 1). We have considered 
the example M = 4, with h = m = uj = 1 and c = 0.5. Then the 8 eigenvalues =bxx are given 
by i^y'p-K, K = 0,1,2,3, with p > 3 and 7 = (5^/2 + 4^/3 - 2)/9 follows from dimi . 
Let us also choose a value for the representation label p: p = M = A. In Figure 3 we give the 
position probability distributions for a number of stationary states w{9)^ namely for 9 = (0, 0, 0, 0), 
(1,0,0,0), (0,1,0,0), (0,0,1,0), (0,0,0,1), (1,1,0,0), (1,1,1,0) and (1,1,1,1). 

Note in these plots that - in agreement with the previous paragraph - only the probabilities 
for ^xk with K = \9\ ox K = \9\ — 1 are nonzero. For 9 = (0,0,0,0), corresponding to the 
highest energy state, the oscillators can be detected only in positions corresponding to the largest 
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eigenvalues xq or —xq. For 6 = (1, 1, 1, 1), corresponding to the lowest energy state, the oscillators 
can be detected only in positions corresponding to the smallest eigenvalues X3 or —X3. The four 
plots with |0| = 1 give some probability distributions in which zizXQ and zizxi are involved. Note that 
for 10(1, 0, 0, 0) and w{0, 0, 1, 0), two stationary states for which the energy level is the same by (|5.2jl 
and (|5.3|1 . also the probability distributions coincide. According to 1)5. 2|) . -E(o,o,o,i) > -^(0,0,1,0) = 
-^(0,1,0,0) > -^(1,0,0,0)- Foi^ the highest of these three energy levels, the probability of detecting the 
oscillator in zizXQ is larger than detecting it in =bxi; for the lowest of these three levels, it is vice 
versa. As |^| increases (thus Eq decreases), the probabilities indicate that the oscillator deviation 
from its equilibrium position also decreases. 

6.4 Coupling of position probability distributions 

In the previous subsection we considered, for a fixed stationary state w{9), the position probabilities 
of the rth oscillator. Due to the symmetry of the system, these probability distributions are 
independent of r. It will be interesting to approach the position probabilities from a different point 
of view. For this purpose, let us assume that the system is in a fixed eigenstate of qi with eigenvalue 
X, say V'i,a:,g- Let us also consider another oscillator r 7^ 1, and the expansion of fpi,x,g in terms of 
the eigenvectors of qr'- 

i^l,x,g = ^M^V'r.s/.h. (6.21) 

y,h 

Then 

EI^m'/ (6-22) 

h 

is the probability of detecting the rth oscillator in the position y (corresponding to the eigenvalue 
y of qr) when the first oscillator is in the state '4'i,x,g- Averaging this out over the multiplicities g 
(if present), thus yields the probability of detecting the rth oscillator in position y when the first 
oscillator is in position x. 

Let us again look at an example of such probability distributions. We shall consider the same 
data as before: M = 4, h = m = u; = l and c = 0.5. First, assume that oscillator 1 is in its 
highest possible position +xq, so the system is in the state ■0i,+xo,i- Then, we can compute the 
probabilities 1)6. 22() . for r = 2, 3, 4 and for y = ±xk (K = 0, 1, 2, 3). These probabilities are plotted 
in Figure 4(a). 

Note that the extreme position of oscillator 1 has a strong influence on the possible positions of 
oscillator 2, a weaker influence on the possible positions of oscillator 3, and again a stronger one on 
those of oscillator 4 (this last is due to the periodic boundary conditions p.2j) . oscillator 4 behaves 
as if it is just to the left of oscillator 1). 

We have also considered a second example, when the first oscillator is in an eigenstate with 
eigenvalue —xi, i.e. the system is in a state V'l.-xi.g- The position probability distributions for the 
other oscillators r = 2,3,4 are plotted in Figure 4(b), where we have averaged out over g (here, 
g = 1,2, 3). Thus in Figure 4(b) an answer is given to the following question: suppose we make a 
measurement of the position of the first oscillator, and detect it in —xi, what are in that case the 
probabilities of finding the other oscillators 2, 3, and 4 in one of their positions ixi^-? 

7 On the spectrum of H and position operators in atypical repre- 
sentations W{p) 

In this section, we very briefly discuss what happens when working with atypical representations 
W{p), i.e. when p G {0, 1, . . . , M — 1}. We concentrate on the spectrum of the Hamiltonian H and 
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on that of the operators and q^. Recah that the representation space in the atypical case is a 
truncation of that in the typical case, discarding basis vectors w{9) for which \9\ > p. Thus one 
expects a close connection between the spectra of the different operators in the typical and atypical 
case. 

It is clear that each basis vector w^O) is still an eigenvector of H, with eigenvalue HEq, with Eg 
given by (|5.2|) . So the spectrum of H in the atypical case is nothing but a truncation of that in 
the typical case where the higher eigenvalues are retained (of course, the actual values are different 
because of the different value for p). As an example, for M = 4 and p = 1 the dimension of the 
representation space is 

and the only four eigenvalues hEg of H are 

-5^(0,0,0,0) = P = Pi + P2 + + Pa, -£'(0,0,0,1) = Pa, 
-5^(1,0,0,0) = Pi = -E'(o,o,i,o) = Ps, and -E(o,i,o,o) = P2, 

where Pi > P3 = Pi > P2 ior < c < cq. These are the four topmost energy levels depicted in 
Figure 2(a). 

In the typical case q"^ has eigenvalues x'j^ = :^^{p — K)'y, with < K < M — 1. From 
equations (|6.1() , ()4.3I) , ()4.7I) and (|4.8() it follows immediately that an eigenvector of q'^ in the typical 
case is also an eigenvector of in the atypical case provided that it is a linear combination of basis 
vectors w{9) with \6\ < p. When K < p one has the same set of eigenvectors as in the typical case, 
arising from both and ■u^ and the multiplicity of x\ is 2(*^^). However, when K = p, Xp = ^ 
and only the vectors arising from Ur remain (the vectors arising from Ur would have \6\ = p + 1 
which is impossible in an atypical representation). The multiplicity of eigenvalue Xp = is thus 
(^-^). Consider the case M = 4 and p = 1; each operator q^ has two eigenvalues namely 0, with 

multiplicity {^\^) = 3 and 4^^(1)7 with multiphcity 2(g) = 2. 

For the position operators q^ finally, it is seen as before that ±rcx with {) < K < p — 1 are 
eigenvalues of qr each with multiplicity (^^^) . Besides these eigenvalues, there is also the eigenvalue 

= with multiplicity {^^p^)- So in the atypical case it is possible to "detect" an oscillator in its 
equilibrium position, in contrast with the typical case. 

It is worth giving some further details for the representation VF(1) (so p = 1), for general 
M-values. This representation has dimension M + 1, with basis vectors u'(O) and w{V) (in the 
notation of ()6.4() ). with j = 1, . . . , M. Each position operator q^- has spectrum {— xq, 0, +xo}, with 

multiplicities {1, M — 1, 1} respectively, where xq = ^/^^- Herein, 7 is given by ()6.11|) : in fact it 
will be useful to introduce the notation 

7fc = ^, fe = l,...,M, (7.1) 

and thus7 = ^^^7f. 

In this case, it is not difficult to construct explicitly a set of orthonormal eigenvectors of qr- In 
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the notation of subsection 6.3, it is given by: 

M 



^^0)-y ^e^u;(F), (7.3) 

1 / V— -^ 27rij-j 1 27rir(o + l) , , \ 

^^.,0,3 = , , E 2^ ' 2 e— ^(1^) e^r^w{l3^') , (7.4) 



^— — ^ + ^ \ 7i + • • • + Tg ■ ■ 79+1 

where g = 1,2, . . . , M - I. 

Now it becomes simple to compute some position probabihties. Fohowing (|6.2UI) . one finds for 
the state w{0) with 6^ = (0) = (0, . . . , 0): 

P((0), r, +xo) = P((0), r, -xq) = ^, P((0), r, 0) = 0. 

Hence in the highest energy state of VF(1), the oscillators can be detected only in the positions ±xo 
and not in 0. Similarly: 

P((P),r,+xo) = P((l^),r,-xo) = ^, P((l^),r,0) = 1-^. 
Also the probabilities (|6.22|) can be computed. One finds, for example: 

^i,+a;o,i - ^ Z^7j COS j, ^i,+xo,i " 7j sm J, 

and thus 

As before, the quantities Sgl^i+xoil^ l^i+xo'il^ describe the probabilities of 

detecting the rth oscillator in the position +xq, or —xq respectively, when the first oscillator is 
in its highest position +xo. 

With the given probabilities, one can consider a final illustration. If the system is in the state 
V'i,+xo,i) the average position of each oscillator r is given by 

-\-Tn\ /l'''+^0,l |2 I n \^ I Ar,0,9 |2 _ i Ar,-xo,l |2 

g 

Using (|7.5j) . this simplifies to 

\ 1 2 /27r(r - l)j 
{(lr)^,,+.„, = ^o-2^7i cos( — ). 

In Figure 5 we plot the average position of each oscillator r in this state. So in this figure one can 
see the effect of having the first oscillator in its highest position +xq on the average position of the 
other oscillators. 
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8 Conclusion and outlook 



We have examined some properties of non-canonical solutions of the quantum system determined 
by the Hamiltonian These solutions arise from a Wigner Quantum Systems approach, where 

the quantization conditions are weaker than the canonical commutation relations, thus allowing 
more types of solutions. 

The solutions studied here are found by identifying certain linear combinations of position and 
momentum operators with generators of the Lie superalgebra 0l(l|Af). We have shown that this is 
always possible, but that the solutions are corresponding to the compact form u(l|M) of 0[(1|M) 
only if the coupling constant c is sufficiently small. 

The physical Hilbert spaces in which the states of the system are described then correspond to 
unitary representations of 0[(1|M). In this paper, we have considered only a simple class of such 
unitary representations, the so-called Fock spaces W{p). This class of representations turns out to 
be already sufficiently resourceful in order to exhibit some fascinating physical properties of the 
solutions. Of special interest is the feature of having only a discrete spectrum for each oscillator 
position operator. At first sight, this is somewhat unusual. On the other hand, our analysis of 
probability distributions for position operators has shown effects that are reminiscent of canonical 
results. 

This paper presents only the first results for this quantum system consisting of a one-dimensional 
chain of coupled harmonic oscillators in the Wigner Quantum Systems approach. There are still 
many open problems or new aspects to be studied. For example, it is clear that the system ()2.29|1 
has also solutions outside g[(l|Af). For instance, if M = 2n is even, then one can construct an 
algebraic solution by means of the direct sum Lie superalgebra 0t(l|2) © • • • ©g[(l|2) (n copies). In 
this case, the unitarity conditions following from the form u(l|2) © • • • ©u(l|2) imply no conditions 
on the coupling constant c, and it would be interesting to study the system from this point of 
view. Furthermore, it would be worth investigating whether p. 29)1 has also solutions related to 
orthosymplectic Lie superalgebras pIT| . 

But even if we restrict for the moment our attention to the g[(l|M) solutions given here in 
Section 01 work remains to be done. In particular, one should also consider other classes of unitary 
g[(l|M) representations [^Hl and investigate the corresponding physical properties. 

Finally, and this a more technical question, the explicit construction of orthonormal eigenvectors 
of the position operators qk is lacking. Although this can be done numerically for any given M and 
a given set of parameters, we have at the moment no closed form expressions for these eigenvectors. 
We hope to find such forms, as they would allow us to draw some general conclusions regarding 
position probability distributions. At this moment, the last conclusions in Section 6 are based upon 
observations of examples rather than upon analytic formulas. We expect to return to some of these 
remaining questions in future publications. 
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Figure 3: Position probability distributions P{6,1,±xk) for a number of ^-values, for M = 4. 
The vertical axis gives the values of P{9, 1, ±xk)- On the horizontal axis one finds the possible 
eigenvalues ±xk with K = 0,1, 2, 3. The values of the parameters are described in the text. 
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Figure 5: Average positions of the M oscillators, M = 10, when the first oscillator is in the position 
+xq, in the representation VF(1). Here, we have taken c = 0.12lj^. The horizontal axis labels the 
M oscillators; the vertical axis gives the average position in units of xq. 
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